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Abstract: 


Similarity  laws  for  the  pressure  spectrum  in  turbulent  flows 
are  developed.  Comparable  to  the  -5/3  power  law  for  the  velocity 
spectrum  in  the  inertial  subrange  there  is  a -7/3  power  law 
relationship  for  the  pressure  spectrum.  Although  the  law  wave- 
number  spectrum  scales  with  relatively  easily  measured  flow  field 
properties,  further  measurements,  especially  of  the  turbulent 
dissipation  are  necessary  to  quantify  the  Holomogorov  constant 
which  appears  in  the  high  wave  number  relationship.  Presently 
the  estimates  for  this  quantity  vary  between  2.25  and  5.6. 
Measurements  of  the  pressure  spectra  in  the  mixing  layer  of  an 
axisymmetric  jet  confirm  the  predicted  scaling  laws  at  high 
and  low  wavenumber . 
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Introduction 

The  idea  of  an  equilibrium  range  in  the  turbulent  energy  spectrum 

was  first  proposed  by  Kolmogorov  (1941  a,b)  who  argued  that  the  parts 

of  the  spectrum  whose  time  scales  were  much  smaller  than  that  for  overall 

decay  might  exist  in  statistical  equilibrium  (c.f.  Batchelor  [1953] ) . 

The  only  characteristic  parameters  for  this  equilibrium  range  are  the 

viscosity,  V,  and  the  rate  of  transfer  of  spectral  energy,  £.  The 

existence  of  such  a spectral  equilibrium  range  was  shown  to  be  exptected 

* 

when  the  turbulence  Reynolds  number,  u £./V,  was  much  greater  than  one; 
in  particular,  (u  £/v)^^  » 1 (c.f.  Batchelor  [1953],  Tennekes  and  Lumley 
[1972]). 

For  sufficiently  large  separation  of  the  energy  containing  range  (the 
range  dominating  the  decay  dynamics)  and  the  dissipation  range  it  was 
further  postulated  that  an  inertial  subrange  should  exist  depending  only 
on  the  rate  of  transfer  of  spectral  energy,  e.  The  existence  of  the 
inertial  subrange  is  contingent  on  even  larger  turbulent  Reynolds  numbers 
than  before;  (r~ )^®  » ^ (c.f.  Tennekes  and  Lumley  [1972],  Batchelor 
[1953]). 

These  ideas  were  applied  with  great  success  to  a variety  of  turbulent 
flow  properties.  The  best  known  result  among  these  is,  of  course,  that 
for  the  turbulent  velocity  spectrum  where  the  form  of  the  spectrum  in  the 
inertial  subrange  was  shown  to  be  given  by 


E(k) 


£2/3  k-5/3 


* 2-.  ,,,  2 2 2. 

u =l/3(u^  + U2  + uy 


(1) 
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where  E(k)  is  the  three-dimensional  energy  spectrum  obtained  by  averaging 
the  directional  spectrum  over  spherical  shells  of  radius  k,  k is  the 
wavenumber,  and  a is  a 'universal'  constant.  As  Kraichnan  (1974)  has 
observed,  the  spectral  predictions  from  universal  equilibrium  theory  have 
had  "an  embarrassment  of  success",  even  in  cases  where  they  would  not  have 
been  expected  to  work.  While  the  universal  equilibrium  arguments  are 
strictly  valid  for  only  homogeneous  flows,  the  use  of  the  local  value  of 
e appears  to  extend  the  validity  of  the  arguments  to  a variety  of 
inhomogeneous  flow  situations  even  when  the  scale  of  the  inhomogeneity 
is  the  same  size  as  the  energy  containing  eddies. 

In  1953,  Batchelor  commented  that  the  pressure  spectrum  "Tr(k)  is 

-7 /3 

readily  found  to  be  Tt(k)^k  for  k » 1/2.".  This  comment  appears  to 
have  been  almost  universally  unnoticed  by  many  recent  investigators  of 
the  turbulent  pressure  field.  Moreover,  the  entire  subject  of  similarity 
scaling  for  the  pressure  spectrum  has  not  been  broached  in  spite  of 
adequate  experimental  data. 

In  this  paper  the  similarity  scaling  laws  for  the  pressure  spectrum 

at  both  high  and  low  wavenumber  are  reviewed.  Following  Tennekes  and 

-7/3 

Lumley  [1972],  the  k range  is  seen  to  result  from  an  asymptotic 
matching  of  the  high  and  low  wavenumber  similarity  spectra  as  the  Reynolds 
number  becomes  large.  Predictions  are  compared  with  recent  experimental 
data  and  the  universal  constant  is  evaluated. 

The  Small  Scale  Spectrum  (large  k) 

For  wavenumbers  large  compared  to  that  characterizing  the  energy 
containing  range  we  expect  the  universal  equilibrium  theory  to  apply. 
Accordingly,  V,  e and  p must  be  the  major  parameters  and  we  have 
Ti=n(k,e,v,p) . The  only  possible  dimensionless  form  is 
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ir(k)  = TT(k) 

2 7/4  3/4  22 

p v e pun 


- Vkn) 


(2) 


3 1 / A 

where  n is  the  Kolmogorov  microscale,  n=(v  /e;  , and  where  u is  the 

1/4 

Kolmogorov  velocity  u=(ve) 

In  view  of  the  local  isotropy  of  small  scale  turbulence  and  the 
success  of  the  universal  equilibrium  scaling  for  the  velocity  spectrum 
we  expect  Equation  (2)  to  be  valid  for  a large  variety  of  turbulent  flows. 

If  the  turbulent  Reynolds  number  is  sufficiently  high,  we  expect  that 
there  exists  an  inertial  subrange;  that  is,  a region  in  which  the  only 
major  parameter  is  e.  Thus  we  have  7r=7r(e,k,p)  and  the  only  possible 
dimensionless  combination  yields 


ir(k)  = a 

P 


2 4/3. -7/3 

p e k 


(3) 


This  was  first  obtained  by  Batchelor  [1953] . 

If  in  the  inertial  subrange  we  associate  the  pressure  fluctuation  at 

a given  wavenumber  with  the  velocity  at  that  wavenumber  we  can  obtain  a 

relation  between  a and  the  corresponding  coefficient  for  the  velocity 

1/3  -1/3 

spectrum.  We  take  u^E  k on  dimensional  grounds  where  u^  is  the 
velocity  of  "eddies"  of  size  1/k.  The  turbulent  energy  spectrum  is  then 
given  by 

2 

Uk 

E(k>  - T 


2/3  , -5/3 
a e k 


(4) 
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It  follows  that 


“k 


- e1/3  k-1/3 


The  pressure  spectrum  is  given  by 


tt  (k)  = 


2 A 
P “l. 


-a  t4/3k-7/3 
P 


from  which  it  follows  that 


(5) 


(6a) 


(6b) 


“P  - “2  • <7> 

For  moderate  Reynolds  numbers  a=1.5  and  therefore  we  expect  0(^=2. 25. 

The  reasoning  above  can  be  extended  to  derive  an  analytical  form 
for  the  entire  wavenumber  spectrum  in  the  universal  equilibrium  range. 

One  of  the  most  useful  forms  of  the  velocity  spectrum  that  has  been 
proposed  is  due  to  Pao  (1965)  and  is  given  by 

E(k)  = a £2^3  k "^3  exp  [-(3/2)  a(kn)^3]  • (8) 

Using  this  for  uk  in  Equation  (A)  and  substituting  into  Equation  (6a)  we 
obtain 


jr(k) 


2 3/A  7/4 
p e v 


a2(kn)  7^3  exp  [-3a(kn)^3] 


(9) 
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This  form  for  the  pressure  spectrum  is  shown  in  Figure  (1) . 

The  deficiencies  of  Pao’s  spectrum  near  k^l/n  are  well  known  and 
are  equally  applicable  here.  Equation  (9)  does  have  the  advantage  that 
it  is  of  closed  form  and  readily  integrable. 

The  Large  Scale  Spectrum 

The  large  scale  spectrum  be  definition  is  dominated  by  those  turbu- 
lent eddies  which  contain  most  of  the  turbulent  energy  and  whose  character- 
istic times  are  the  same  order  as  the  decay  time.  Since  the  process  by 
which  energy  is  added  to  the  turbulence  is  geometry  dependent  we  should 
not  expect  a universal  scaling  as  in  the  equilibrium  range  but  at  best 

similarity  spectra  which  are  valid  only  for  a particular  class  of  flows. 

~2 

On  dimensional  grounds  p , the  local  mean  square  turbulent  pressure 
fluctuation,  and  £,  a characteristic  length  for  the  turbulence,  might 
provide  a useful  scaling.  We  have 


= F (k£) 
—y  PP 
A*  P 


(10) 


The  characteristic  length,  i,  remains  to  be  specified.  Since  we  are 

seeking  a length  characteristic  of  the  decay  process,  we  will  choose  the 

3 

length  defined  by  e=u  /l  where  e is  the  local  rate  of  decay  of  turbulent 
energy  and  u is  the  local  rms  turbulent  velocity. 

In  an  incompressible  flow  the  pressure  is  given  by 


•%2  , , 

2 8 Vi 

V p(x,t)  = -p  1 


9x^3xj 


(ID 
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where  is  the  total  velocity  vector.  This  can  be  solved  to  yield 

<\2__ t 

= fa 

space 

It  is  clear  that  the  pressure  depends  on  not  only  the  flow  at  the  point 
in  question  but  on  the  flow  at  a distance.  Thus,  there  is  no  reason, 
a priori,  we  should  expect  that  the  length  Z will  accomplish  the  scaling 
we  seek.  It  might  be  noted  here  that  the  same  considerations  are  valid 
for  the  turbulent  energy  spectrum  where  e must  be  reasonably  constant 
over  distance  of  at  least  Z if  the  large  scale  turbulence  structure  is  to 
be  scaled. 

In  general,  more  information  is  known  about  the  turbulent  velocity 

field  than  about  the  turbulent  pressure  field.  It  would  be  very  useful 

if  the  pressure  spectra  could  be  scaled  in  terms  of  velocity  parameters 
~2 

alone  instead  of  p . We  propose  similarity  relations  for  the  large  scale 
spectrum  based  on  p,  u,  and  Z and  attempt  to  indicate  where  the  relations 
might  be  useful.  From  dimensional  considerations  we  can  write  immediately 
for  the  large  scale  spectrum 


o u . u . 

1 J 


3y,3y_.  I y-x  I 


d3y 


(12) 


TT(k) 

2 4 , 

p u 'i 


= F (k £) 
PP 


(13) 


To  examine  the  range  of  validity  of  Equation  (13)  we  will  look  at  the 
dependence  of  the  turbulent  pressure  field  on  the  turbulent  velocity  field. 
For  a homogeneous,  isotropic  turbulence  Batchelor  [1933]  was  able 


to  show  that 
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f 4 

ir(k)  = E(k')E(|k-k’|)  -5-1-  °T  dk'  . (14) 

8tt  ‘ ~ ~ | k-k ' | 

This  relation  was  based  on  an  assumption  of  normality  for  fourth  order 
moments  and  is  not  expected  to  be  valid  for  high  wavenumbers.  From 

2 

Equation  (14)  Batchelor  argued  that  the  principal  contribution  to  p came 

from  wavenumbers  near  the  maximum  of  E(k) , the  velocity  spectrum.  Since 

the  maximum  of  E(k)  is  at  about  1 /l  we  should  expect  Equation  (13)  to  work 

for  this  restricted  class  of  turbulent  flows. 

Kraichnan  [1956]  treated  a slightly  more  general  class  of  flows  than 

Batchelor  and  included  the  effects  of  anisotropy  and  a uniform  shear. 

The  effect  of  anisotropy  was  found  to  weaken  slightly  the  mean  square 

pressure  fluctuation  for  a given  turbulent  energy.  For  the  case  of  an 

isotropic,  homogeneous  turbulence  acted  upon  by  a uniform  and  steady  mean 

shear,  Kraichnan  showed  that  there  were  two  contributions  to  the  mean 

~2 

square  pressure,  a contribution  from  the  turbulence  alone,  p T_T>  and  a 

~2 

contribution  from  the  turbulence-shear  interaction,  p Under  the 

assumption  of  normality  for  fourth  order  moments  and  an  exponentially 
decaying  spatial  correlation  for  the  velocity  field,  these  terms  were 
found  to  be 


P2t_x  ^ U/2puV 


and 


(15) 


2 „ 2 2„2  ,du, 2 

'l-s'1'1  W 


Arndt,  et  al.  [1974]  and  others  have  observed  reasonably  good 
agreement  with  these  predictions  in  jet  mixing  layers.  The  predictions 
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of  Equation  (15),  however,  do  not  seem  adequate  in  the  adjustment  region 
of  jets  or  in  more  general  classes  of  flows  such  as  a body-perturbed  jet. 

As  noted  by  Kraichnan,  isotropic  homogeneous  turbulence  subjected  to 
steady  uniform  mean  shear  is  a highly  artificial  example  since  such  a 
flow  would  quickly  evolve.'  The  success  of  Kraichnan' s theory  in  mixing 
layers  is  probably  due  to  the  fact  that  the  turbulence  is  always  in 
local  scale  equilibrium;  that  is  the  ratio  (&/u) (du/dy)  - constant.  If 
this  is  the  reason,  then  Kraichnan's  result  is  forced  by  dimensional 
considerations  and  probably  will  not  be  valid  in  flows  involving  more  than 
one  time  scale  as  in  the  wall  region  of  a boundary  layer  or  in  a perturbed 
jet.  This  appears  to  be  the  case.  It  is  interesting  to  speculate  whether 
Kraichnan's  results  could  be  extended  to  a more  general  class  of  flows  by 
using  Lumley's  [1967a]  concepts  of  limited  spatial  awareness  and  limited 
temporal  memory. 

In  flows  which  are  characterized  by  a single  time  and  length  scale 
we  expect  Equation  (13)  to  provide  scaling.  For  more  general  flows,  it 
will  be  necessary  to  revert  to  the  form  proposed  in  Equation  (10)  which 
requires  only  that  Z be  the  appropriate  length. 

The  Turbulence-Turbulence  Interaction  Spectrum 

To  date  no  satisfactory  method  for  experimentally  separating  the 
pressure  fluctuations  due  to  the  turbulence-turbulence  interaction  and  the 
turbulence-mean  shear  interaction  has  been  proposed.  In  the  paragraphs 
below  a spectral  form  for  the  turbulence-turbulence  pressure  spectra  is 
proposed.  It  is  suggested  that  spectral  deviations  from  this  form 
represent  the  turbulence-mean  shear  interaction. 


H 
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By  assuming  that  the  work  done  by  the  mean  strain  rate  per  unit 
wavenumber  and  per  unit  time  is  determined  only  by  local  quantities, 
Tennekes  and  Lumley  [1972]  obtained  a low  wavenumber  velocity  spectrum  for 
turbulent  flow: 


E(k)  = a e2/3  k 5/3  exp  [-3/2ir  $ ct1/2(k£)  4/3] 


(16) 


where  a is  the  Kolmogorov  constant  and  £ is  determined  by  £=u  /£.  The 

2 

spectrum  was  required  to  integrate  to  the  kinetic  energy  3/2  u which 
yielded  3=0.3.  The  model  was  not  expected  to  be  valid  for  k£  <<  1. 

Since  the  hypothesis  of  dependence  on  local  parameters  has  already 
been  made,  we  can  apply  arguments  similar  to  those  preceding  Equation  (9) 
and  obtain  for  the  low  wavenumber  spectrum 


Fpp(k£)  = a2  (k£)-7/3  exp  [-3tt  3 a1/2  (k£)  4/3] 


(17) 


or 


tt (k)  = a2p2eA/3  (k£)_7/3  exp  [-3tt  3 a1/2  (k£)_4/31 


This  is  shown  graphically  in  Figure  (2) . 

The  Inertial  Subrange 

The  dependence  of  the  existence  of  an  inertial  subrange  on  turbulent 
Reynolds  number  can  best  be  illustrated  by  a technique  used  by  Tennekes 
and  Lumley  [1972]  for  the  velocity  spectrum.  We  have  proposed  similarity 
laws  for  the  spectrum  at  both  high  and  low  wavenumber.  As  the  turbulent 
Reynolds  number  becomes  large,  the  energy  containing  range  and  the  dissi- 
pation range  become  more  widely  separated.  Therefore,  we  ask  whether  there 


\ 
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exists  a region  of  common  validity  for  both  spectral  forms  if  we  take 

the  limits  R-**5,  k£-x»,  and  kry*0  simultaneously. 

-3/4 

It  is  well  known  that  qJ/^R  (Tennekes  and  Lumley  [1972]).  To 
take  the  limits  we  let  kA^Rn  where  n>0.  In  order  for  kq->-0  as  R-*00 
we  must  have  k =ki.(n/£)^Rn  3^-*0.  Therefore,  any  n satisfying  0<n<3/4 
will  suffice. 

Substituting  these  relations  into  Equations  (2)  and  (13)  we  obtain 


pV'V/4  £ <r"-3/4>  - pV*  F (R") 

PP  PP 


which,  using  £=u  /l,  becomes 


f (Rn"3/S  = R7M  F (Rn)  . 
PP  PP 


(18) 


(19) 


This  must  be  satisfied  for  any  n in  the  interval  0<n<3/4.  The  solution 
is  readily  obtained  as 


f (kn)  = a (kn) 
PP  P 

F (kq)  = a (k£) 
PP  P 


-7/3 

-7/3 


(20) 


Thus  there  is  a region  of  overlap  and  it  corresponds  to  the  equilibrium 
range  found  previously. 

The  Mean  Square  Pressure  Fluctuation 

The  mean  square  pressure  fluctuation  is  obtained  from  the  spectrum 


by 
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— r 2 2 

2 p u 


u • £ 
£ 


F (k£)  dk£ 
PP 


For,  geometrically  similar  flows  the  integral  is  nearly  a constant  and 
we  can  write 


p = A1  • p u 


The  slight  dependence  of  A^  on  Reynolds  number  results  from  the  increasing 
-7/3 

extent  of  the  (k£)  range  with  Reynolds  number  when  the  low  wavenumber 

scaling  is  applied. 

A value  for  A^  can  be  obtained  for  high  Reynolds  numbers  by  integrating 
the  spectral  form  of  Equation  (17) . 


w 

?2  = p2e4/V/3a2  (k  )'7/3  exp  [-3-IT  B a1/2  (k£)_4/3]  dk£ 


£!1 

4T7&  £ % P 


Using  a=1.5  and  £=u  we  obtain 


2 ~ n /a7  2 22 
p s=  0 .487 p u 


Hinze  [1959]  calculated  p'  from  an  assumed  form  of  the  velocity 

correlation  f(r)  = u-^(x)  u-^x+r)  for  isotropic  turbulence.  Using  simple 

~2  2 ~22  , 

exponential  decay  (high  Reynolds  number)  he  obtained  p - 0.49p  u and 

~ 2 2~ 2^ 

by  using  a Gaussian  correlation  (low  Reynolds  number)  p - l.Op  u 
Batchelor  [1951]  used  the  velocity  correlation  measured  in  a grid  turbulence 


J 
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2 2 

to  obtain  p - 0.34  p u . The  estimate  obtained  here  and  the  high 
Reynolds  estimate  of  Hinze  are  in  excellent  agreement  with  the  results 
of  Uberoi  [1954]  who  obtained  ^'-0.7  pu'  . This  agreement  with  Uberoi's 
results  is  particularly  gratifying  for  three  reasons:  first,  there  are 

no  adjustable  constants  in  this  theory;  second,  Uberoi's  results  were 
obtained  in  a nearly  isotropic  turbulence  where  no  turbulence  shear 
interaction  term  exists;  and  third,  the  results  were  inferred  from 
velocity  measurements,  thus  avoiding  the  question  as  to  whether  the  probe 
was  measuring  pressure  fluctuations. 

The  Mean  Square  Fluctuating  Pressure  Gradient 

The  mean  square  fluctuating  pressure  gradient  is  given  by 

00 

(Vp)2  = | k2  7T(k)  dk  . (26) 

o 


The  integrand  will  rise  rapidly  for  low  wavenumbers  and  peak  near  the 

low  wavenumber  end  of  the  inertial  subrange  beyond  which  it  will  roll 

-1/3  -1/3 

off  as  k . Thus  most  of  the  contribution  will  come  from  the  k 

range  for  high  Reynolds  number  flow.  It  is  interesting  to  note  that 

this  spectrum  does  not  peak  at  high  wavenumbers  as  does  the  velocity 

derivative  spectrum. 

To  evaluate  the  mean  square  pressure  gradient  for  large  Reynolds 
number  we  use  the  high  wavenumber  spectrum  of  Equation  (9)  and  obtain 


,n  ,2  2,  3,  .1/2  2 

(Vp)  = p (e  /v)  a 


(kn)  exp  [-3a  (kri)^^]  dkn 


(26) 


o 
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2/3 

Making  the  variable  change  y = (kn)  this  can  readily  be  integrated  to 
yield 


(Vp)2 


/3tT  3/2  2 
—r~  a p 


(e3/v) 


1/2 


(27a) 


or  using  01=1.5  we  obtain 


(Vp)2  » 4.0  p2  (e3/v)1/2 


(27b) 


Batchelor  [1951]  obtained  a similar  result  from  an  approximate 
spectral  model. 

The  One  Dimensional  Spectrum 

While  the  three-dimensional  spectrum  TT(k)  is  useful  for  dynamical 

considerations,  it  is  the  one-dimensional  spectrum  F3  (k.)  that  is 

PP  1 

usually  measured.  Their  relationship  is  described  below. 

The  pressure  covariance  for  homogeneous  turbulent  flow  is  defined 
by 


Rpp(~)  = P^~  + 


(28) 


and  its  Fourier  transform,  the  pressure  spectrum  is 


* (ki»k9»k-,)  = n 

PP  1 2 3 (27t)3 


00 
f f 


„ . . -ik*r  , 

R (r)  e ~ ~ dr 
PP  ~ 


(29) 


The  one-dimensional  pressure  spectrum  F^(k^)  is  defined  by 
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F*  (k  ) 
PP  1 


%p^kl,k2’k3^  dk2dk3 


(30) 


The  three  dimensional  spectrum  ff(k)  is  obtained  from  ^^(k)  by  averaging 
over  spherical  shells  of  radius  k;  that  is 


1T(k)  = 

47Tk 


$ (k)  do  (k) 

PP  - 


(31) 


where  do(k)  is  an  element  of  area  on  a sphere  of  radius  |k|  = k. 

If  the  turbulence  is  isotropic  there  exists  a relation  between 
ir(k)  and  F^  (k,)  (c.f.,  Hinze  [1959]  and  Lumley  [1970]): 

pp  1 


or 


fJLoo 

PP  1 


= 1/2 


1(10 


dk 


TT(k)  = -2k 


3Fin(k) 

E£ 

9k 


(32a) 


(32b) 


It  is  clear  from  Equation  (32b)  that  a power  law  region  of  TT(k)  will 

also  appear  as  the  same  power  law  in  F^  (k.  ) . Thus  the  one-dimensional 

PP  1 
-7 /3 

velocity  spectrum  will  also  exhibit  a k region.  We  obtain  for  the 


inertial  subrange 


f1  (k.)  = 3/14  a k,  7/3  . (33) 

pp  1 pi 

The  one-dimensional  pressure  spectrum  analog  of  Pao's  velocity 
spectrum  can  be  obtained  by  substituting  Equation  (9)  into  Equation  (32a) 
and  integrating.  In  dimensionless  form 
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■ 1/2  «2 


r-io/3  , . r-4/3,  ,r 
k exp  [-3  ak  ] dk 


where 


(34) 


_PJL 


PP 


2 3/4  7/4 
p e v 


k 


1 


k = kn 


(35) 


Following  Pao  [1965]  we  introduce  the  variable  £ = k^/k  and  obtain 


<*i>  ' 1/2  °2  i'1‘7/3 


C4/3  exp  [-3a  (^z^)473]  dC  • (36) 


This  has  been  integrated  numerically;  the  results  are  tabulated  in  Table  T 
and  shown  graphically  in  Figure  (3) . 

A similar  exercise  can  be  carried  out  for  the  large  scale  spectrum 
with  the  following  integral  resulting: 


FpP<u> 

p2u,4Jl 


1/2  a2  (k1JO~7/3 


£4/3  exp[-3TrBa1/2 


k.  -4/3 

i-f)  ] dC 


(37) 


This  has  also  been  integrated  numerically;  the  results  are  tabulated  in 
Table  II  and  shown  graphically  in  Figure  (4) . 

It  should  be  remembered  that  both  of  these  results  assumed  isotropy 
and  therefore  their  applicability  to  experimental  data  is  limited  to  such 
flows.  On  the  other  hand,  in  view  of  the  local  isotropy  of  most  high 
Reynolds  number  turbulent  shear  flows  at  smal)  scales,  the  high  wavenumber 
spectrum  should  find  wide  application. 
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The  Experimental  Results 

An  attempt  was  made  to  collapse  the  considerable  body  of  data  that 
has  appeared  in  the  literature  in  recent  years  (c.f.,  Planchon  [1974], 

Fuchs  [1972],  Arndt,  Tran,  and  Barefoot  [1974]).  Although  all  of  these 
pressure  spectral  measurements  indicated  a -7/3  region,  dimensionless 
scaling  was  not  possible  because  estimates  of  u,  £,  and  e were  not 
available.  Therefore  an  investigation  was  conducted  in  the  mixing  layer 
of  a 12  inch  open  jet. 

The  unsteady  pressure  probe  used  in  this  investigation  was  developed 
by  Arndt  and  Nilsen  [1969].  The  sensitive  element  is  a Bruel  and  Kjaer 
1/8  inch  condenser  microphone  which  is  connected  to  a cathode  follower 
and  powered  by  a B & K Type  2801  Microphone  Power  Supply  Unit.  The  probe 
is  a standard  pitot  static  tube,  0.12  inches  outside  diameter  and  2.5  inches 
in  length.  Four  static  pressure  holes  are  spaced  90  degrees  apart  at  a 
distance  of  5 tube  diameters  from  the  tip  of  the  probe;  the  other  end  of 
the  probe  is  terminated  by  the  microphone.  Damping  material  was  placed 
in  the  tube  to  diminish  a response  peak  at  tube  resonance  frequency. 

The  probe  was  calibrated  by  placing  it  and  a reference  microphone  in 
a reverberent  chamber.  The  calibration  and  a sketch  of  the  probe  are 
shown  in  Figure  (5).  Spectra  were  measured  using  a Federal  Scientific 
Real-Time  Analyzer.  The  effective  bandwidth  was  10  Hz  and  the  effective 
bandwidth-averaging  time  product  was  400.  Spectra  were  plotted  on  a x-y 
plotter  as  log  spectrum  versus  frequency  and  the  data  points  shown  were 
read  from  these  plots.  This  procedure  introduced  errors  estimated  at  20%. 
The  measured  spectra  were  corrected  for  the  temporal  response  of  the  probe 
up  to  a maximum  correction  of  2dB.  No  attempt  was  made  to  correct  for  the 
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limited  spatial  response  which  was  judged  to  be  considerable  at  frequencies 

above  f=u  /10d  where  u is  the  convection  velocity  and  5d  is  the  distance 
c c 

from  the  probe  tip  to  the  static  holes. 

Arndt,  Tran,  and  Barefoot  [1974]  using  the  same  probes  report  that  at 
the  center  of  the  mixing  layer 

p’  = 0.0625  (l/2p  u"2) 

which  corresponds  to 

p'  = 1.26  pu2 


Hence  for  this  flow  it  can  be  assumed  that  the  shear  interaction  term  is 
of  the  same  order  as  the  turbulence  term.  Thus,  this  experimental  situation 
provides  a good  test  for  the  scaling  arguments  of  the  previous  sections. 

Spectral  measurements  were  made  in  the  center  of  the  mixing  layer  at 
locations  of  x/D  = 1,  1.5,  2,  3 and  4 where  D is  the  jet  diameter.  The 
exit  velocities  used  were  19.8  m/s  and  30.5  m/s.  Since  the  exit  diameter 
was  0.305  m,  the  exit  Reynolds  numbers  were  4.0  x 10^  and  6.2  x 10~* 
respectively.  The  previous  work  of  Arndt  et  al.  [1974]  indicated  that  the 
lower  velocity  was  the  smallest  that  could  be  used  without  appreciable 
contamination  from  background  noise. 

There  have  been  several  exhaustive  series  of  measurements  of  velocity 
data  in  the  axisymmetric  jet  mixing  layer  (c.f.  Davies,  et  al . [1963], 
Bradshaw  et  al.  [1964]).  Unfortunately,  the  rate  of  dissipation  of 
turbulent  energy  per  unit  volume,  e,  was  not  measured.  Because  of  the 
extremely  small  probes  that  would  have  been  required  to  accurately 
determine  the  necessary  velocity  derivatives,  these  measurements  were  not 
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carried  out.  It  is  possible,  however,  to  state  the  approximate  functional 
dependence  of  u,  e and  l on  the  distance  from  the  jet  lip,  x,  and  the 
exit  velocity,  U^.  These  are  readily  obtained  as 


from  which  it  follows  that 


SL  ^ x 
U%UE 


Thus  the  data  obtained  in  this  experiment  are  plotted  with  two 
nondimensionalizations  corresponding  to  the  high  and  low  wavenumber 
forms  discussed  earlier: 


and 


f U„x 

n,  P-P—  (-L-) 

23/47/4  2 4 _ v v ; 


_EP 


p e v 


P UE  x 


i „ , , v . 3/4 

kn  ^ kx  (— — ) 

UE 


PP 


2 4 . 2..  4 

P u SL  p U£  x 


kS,  ^ kx 


high  wavenumber 


low  wavenumber 


The  Kolmogorov-like  spectra  (high  wavenumber  scaling)  are  shown  in 
Figure  6.  Wavenumbers  were  calculated  by  using  a fixed  convection  velocity 
equal  to  the  mean  velocity  at  the  mixing  layer  center  (U  = 0.6  U ).  In 
view  of  the  accuracy  of  the  results,  application  of  Lumley's  [1967b] 
correction  to  Taylor's  hypothesis  hardly  seems  justifiable  (less  than  20%). 


JL 


\ 
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In  particular,  the  measurement  of  the  mean  velocity  is  accurate  only  to 
within  5%.  Moreover,  since  measurements  were  taken  at  the  center  of  the 
mixing  layer  where  the  mean  velocity  gradient  is  a maximum,  additional 
error  could  result  from  a slight  mislocation  of  the  probe.  Since  velocity 
enters  to  the  fourth  power',  the  observed  20%  scatter  is  not  unreasonable. 
There  is  some  evidence  of  such  systematic  error  in  the  normalized  plots. 

The  roll-off  at  high  wavenumbers  is  attributed  to  the  spatial 
averaging  of  the  probe.  The  slight  positive  deviation  from  the  -7/3 
law  deviation  is  probably  due  to  background  noise  since  it  progressively 
disappears  as  the  absolute  spectral  level  increases.  The  line  drawn 
through  the  data  is  chosen  to  fit  the  highest  absolute  data  for  which  the 
background  effect  would  be  the  least  (x/D  = 1,  U = 30.5  m/s). 

The  low  wavenumber  plots  of  measured  spectra  are  shown  in  Figure  (7) . 

It  is  clear  that  the  inertial  subrange  also  collapses  on  this  plot.  The 
low  wavenumber  data  do  not  collapse  and  the  spectral  value  is  seen  to 
decrease  with  distance  from  the  origin.  Certainly  Taylor's  hypothesis  is 
not  valid  for  wavenumbers  in  the  range  k£  < 1 (c.f.  Lumley  [1967b]).  In 
particular,  the  eddy  turnover  time  or  life-time  is  of  the  same  order  as  the 
time  required  to  pass  the  probe.  Such  behavior  would  also  be  observed  if 
the  pressure  fluctuations  resulting  from  the  shear  interaction  were  dominated 
by  the  large  eddy  structure  (c.f.  Arndt  and  George  [1974]).  Therefore  this 
disparity  at  low  wave  numbers  is  not  felt  to  contradict  the  low  wavenumber 
scaling  hypothesis. 

Wygnanski  and  Fiedler  [1970]  measured  the  Taylor  microscales  in  the 
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local  isotropy  (which  was  not  confirmed  in  the  work  cited)  one  obtains 

UE 

e = 0.046  — . 

x 

Using  u=0.16  Ug  from  the  same  work  and  confirmed  in  this  study,  one 
obtains  £=0.08x.  When  these  values  are  substituted  into  the  scaling  laws 
and  compared  to  the  measured  pressure  spectra,  the  value  of  Op  is  obtained 
as  5.6.  If  on  the  other  hand  £ is  selected  as  5/v0.04x  the  value  suggested 
earlier  can  be  obtained;  i.e.,  a = 2.25.  In  view  of  the  difficulties 
associated  with  accurate  dissipation  measurements  in  such  flow,  the  latter 
values  does  not  seem  unreasonable. 

Figure  (8)  shows  a plot  of  the  pressure  spectra  normalized  with  exit 
parameters  U^,  D and  Strouhal  number.  It  is  clear  that  the  pressure 
spectra  maxima  do  not  collapse  in  this  type  of  plot  and  the  spectral 
peak  clearly  shifts  to  lower  frequencies  with  increasing  x as  would  be 
expected . 

In  all  plots  it  is  clear  that  for  the  lowest  frequencies  the  spectrum 
increases  quadratically  with  increasing  frequencies.  The  peak  and  slope 
is  taken  to  indicate  that  the  pressure  space  correlation  is  negative 
somewhere  (c.f.  Luraley  [1970]).  This  interpretation  is  consistent  with 
the  observations  of  Planchon  and  Jones  [1974] . 

The  fact  that  the  spectra  display  a limited  -7/3  range  is  interpreted 
here  as  confirmation  that  the  pressure  probe  is  indeed  measuring  static 
pressure.  If  the  probe  were  responding  to  stagnation  pressure,  the  dominant 
fluctuating  term  would  be  the  cross  product  of  the  mean  and  fluctuating 
velocity;  thus  a -5/3  spectrum  would  have  been  observed.  This  is  perhaps 
the  explanation  for  the  data  of  Elliot,  et  al.  [1972]. 
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Summary  and  Conclusions 

Scaling  laws  have  been  proposed  for  both  high  and  low  wavenumber 

pressure  spectra.  Exact  spectral  forms  for  both  high  and  low  wavenumbers 

-7/3 

have  been  derived.  An  inertial  subrange  behaving  as  Tr(k)Mc  has  been 
shown  to  exist.  Further  measurements,  especially  of  the  turbulent 
dissipation,  are  necessary  to  evaluate  the  Kolmogorov  constant  although 
it  appears  to  be  between  the  value  of  2.25  suggested  and  5.6.  Meaurements 
of  pressure  spectra  in  the  mixing  layer  of  an  axisymmetric  jet  were  shown 
to  follow  the  predicted  scaling  laws  at  both  high  and  low  wavenumbers. 
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Figure  7 - Pressure  Spectra  Scaled  with  Large  Scale-Like  Variables  (q,®ho0 
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UE 


Figure  8 - Pressure  Spectra  Scaled  with  Exit  Parameters,  U£,  D 
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